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ABSTRACT 

We derive hydrostatic, radiative equilibrium models for passive disks 
surrounding T Tauri stars. Each disk is encased by an optically thin layer of 
superheated dust grains. This layer re-emits directly to space about half the 
stellar energy it absorbs. The other half is emitted inward and regulates the 
interior temperature of the disk. The heated disk flares. As a consequence, it 
absorbs more stellar radiation, especially at large radii, than a flat disk would. 
The portion of the SED contributed by the disk is fairly flat throughout the 
thermal infrared. At hxed frequency, the contribution from the surface layer 
exceeds that from the interior by about a factor 3 and is emitted at more than 
an order of magnitude greater radius. Spectral features from dust grains in the 
superheated layer appear in emission if the disk is viewed nearly face-on. 


Subject headings: circumstellar matter — radiative transfer — stars: pre-main 
sequence — infrared: stars — accretion, accretion disks — stars: individual 
(CM Aur, GG Tan) 


1. INTRODUCTION 

Excess infrared (IR) emission from T Tauri stars is thought to arise from circumstellar 
disks (Mendoza 1968; Shu, Adams & Lizano 1987; Adams, Lada, & Shu 1987; and references 
therein). The standard disk is modeled as a flat blackbody. In the simplest scenario it 
lacks intrinsic luminosity and passively reradiates the energy it absorbs from the central 
star. This blackbody model yields an infrared spectral energy distribution (SED) of the 
form uFt, oc u"', with spectral index n = 4/3 (Adams, Lada, & Shu 1987; and references 
therein). By coincidence, the steady-state accretional luminosity from an opaque, thin disk 
yields an SED with identical spectral index (Lynden-Bell & Pringle 1974). SEDs measured 
between 3 and 100 pm are usually well htted by power laws with 0 < n < 4/3. However, 
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in contrast to the prediction of the standard model, most sources exhibit flattish spectra 
having n < 3/4 (see, e.g., Rydgren et ah 1984; Rucinski 1985; Rydgren & Zak 1987; Strom 
et ah 1989; Beckwith et ah 1990). 

The failure of the standard disk model spawned a number of alternative proposals. Of 
greatest relevance here is the investigation by Kenyon & Hartmann (1987) of a blackbody 
disk whose surface flares outward with increasing radius as a consequence of vertical 
hydrostatic equilibrium. Flared disks intercept more stellar radiation than flat ones, 
especially at large distances from the star. Other models invoke either an “active” disk 
having a high intrinsic luminosity, or a dusty component in addition to the disk. The 
dissipation of a one-armed spiral density wave as an intrinsic disk heating mechanism falls 
in the hrst category (Adams, Ruden, & Shu 1989; Shu et ah 1990; Ostriker, Shu, & Adams 
1992). The second category includes the proposal that a tenuous dusty envelope surrounds 
the star/disk system and either scatters stellar radiation back onto the disk (Natta 1993), 
or absorbs and re-emits it in the infrared (Calvet et al. 1994). 

We calculate the SED for the passive, reprocessing disk in a self-consistent fashion. 
The basic features and results of our investigation are set forth in §2. Extensions and 
rehnements are considered in §3. Lastly, in §4, we discuss our results in the context of 
observations and summarize unresolved issues. 


2. THEORY 
2.1. Model Assumptions 

We consider a T Tauri star surrounded by a passive disk. The star is modeled 
as a spherical black body of temperature = 4000 K, mass = O.SMq, and radius 
R* = 2.5Rq (see, e.g.. Table II in Beckwith et al. 1990). Our hducial disk has surface 
mass density similar to that of the minimum-mass solar nebula, S = a^u^^So, where oau 
is the disk radius measured in AU and Sq = lO^gcm”^ (Weidenschilling 1977). Dust, 
which is uniformly mixed with the gas, comprises about one percent of the total mass. 
Dust grains dominate the continuum opacity from visible through millimeter wavelengths. 
We take the grains to be spheres with radius r = O.l/xm, mass density pd = 2gcm“^, and 
negligible albedo. Their emissivity, which is unity for A < 2nr\ decreases as Sa ~ (27rr/A)^ 
for A > 27rr. We denote by e ~ {87irkT/hc)^ ~ {T/T^Y fhe average dust emissivity at 
temperature T. For our hducial disk, {3 = 1. The dust opacity at visual wavelengths is 
Kv ~ 400cm^g“^, which implies an optical depth ~ 4 x 10® hducial disk. 

SEDs are computed for disks viewed pole-on. We choose an inner cutoh radius. 
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tti ~ ~ 0.07AU, to mark the condensation boundary of common silicates. The outer 

cutoff radius for flat disks, Oq ~ 2.3 x ~ 2.7 x 10^AU, is fixed to facilitate comparisons 

among different disk models; Oq is comparable to the size of the largest disks seen in 
silhouette against the Orion nebula (McCaughrean & O’Dell 1996). 

Our model assumptions are chosen to make the telling of our story as simple and direct 
as possible. Issues regarding these choices and possible alternatives are dealt with in §3. 


2.2. The Blackbody Disk 

To begin, we review standard relations for the temperature and SED of a blackbody 
disk. The flux of stellar radiation incident upon the disk is (a/2) (i?*/a)^cr T/ for a S> R*, 
where a is the grazing angle at which the starlight strikes the disk. Equating the emitted 
and absorbed fluxes yields the disk temperature^ 
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T 

-L * • 
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The SED is computed as 

Lu = = Stt^z/ f daaB^{Te ), (2) 

Jai 

where B^{T) is the Planck function, and d is the distance to the source. A scaling relation 
for Li, at wavelengths between those that characterize the disk at and Oo is derived as 
follows. Equation (j^) is approximated as S'K'^a^uB^{Te), pByiTe) is replaced by crTg/vr, a is 
related to Tg by equation (P, and Tg is expressed in terms of v by 3/cTg ~ hv. These steps 
yield 


Ly ~ ~ aL*. (3) 

In other words, the fraction of the stellar luminosity L* that is reprocessed to frequencies 
in an octave centered on u is approximately equal to the grazing angle a at the location 
where SkT^ = hv. 


^The subscript e is used to denote effective temperature. 
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2.2.1. Flat Geometry 

A flat disk is one whose aspect ratio (opening angle) is independent of a; we assume it to 
be much less than unity. The grazing angle appropriate to this geometry is a ~ 0.4i?*/a -C 1 
for a/i?* 3> 1 (Kusaka, Nakano, & Hayashi 1970, Ruden & Pollack 1991). Thus Tg from 
equation (|I|) takes the form 


2 /R \3/4 

^ T* 

.OTT/ v a 


( 4 ) 


Application of the scaling relation given by equation (| 
Ly ~ 0.01(i//10^^ Hz)^/^L*. Figure |l| conhrms that Ly 
30 ^ A ^ 1000 /im. 


J) to the flat disk gives 
obeys this relation for 


2.2.2. Hydrostatic Equilibrium 

We retain the assumption that the disk radiates as a blackbody and consider the 
consequences of vertical hydrostatic equilibrium in a gravitational held g = VFz. This is the 
case investigated by Kenyon & Hartmann (1987). 

The disk temperature Tg is still given by equation (Q), but now the grazing angle a 
takes the more general form 


0.47?* d 

a ~-h a— 

a da 
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where H is the height of the visible photosphere above the disk midplane. 

Taking the gas to be isothermal at temperature Tg yields a Gaussian vertical density 
prohle 


where 


n 

— = exp 
no 



( 6 ) 
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( 7 ) 
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Fig. 1.— SED for the flat blackbody disk, with contributions from star and disk identihed. 
The n = 4/3 law is evident between 30 /rm and 1 mm. The turnover near 1 mm is due to our 
truncation of the disk at Oq ~ 270 AU. 
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The temperature Tc is a measure of the gravitational potential at the surface of the central 
star; 


T = 

c — 




kR^ 


8 X 10® K, 


( 8 ) 


with Hg the mean molecular weight of the gas. Under the assumption that the dust to gas 
ratio is uniform throughout the disk, we hud 


H 

~h 



( 9 ) 


where nph is the number density in the photosphere. 

Armed with equations (H)-®, we are ready to evaluate H/a. In the limit of large 
radius where a is dominated by the flaring term, 


H 

a 



4/7 



0.18a 


2/7 
AU • 


( 10 ) 


In writing equation (pUD we have set H/h = A] in reality, this factor declines from about 5 
at oau = 3 to 4 at oau = 10^. 


It follows from equations and (pT|) that a ~ O.OOSa^u + 0.05aAu- Thus a is minimal 
at the transition radius atr ~ 0.4AU with amin ~ 0.05. Beyond the transition radius the 
disk flares until H k, a bX an ^ 270 AU. 


The SED for the flared blackbody disk truncated at as computed from equation 
(y), is shown in Figure |^. At wavelengths between 10 and 100 pm, the disk emission follows 
the scaling relation [cf. eq. (^] Ly ~ cxL* ~ 0.1(z//10^^ Longward of 300 pm in 

the Rayleigh-Jeans regime, the SED varies as 


2.3. The Radiative Equilibrium Disk 

Here we drop the blackbody assumption and determine the SED by application of the 
techniques of radiative transfer. The dehnition of the effective temperature, Tg, remains 
as given by equation (|^. A simplihed description of the disk distinguishes two regions. 
The surface layer contains grains which are directly exposed to light from the central star. 
Variables evaluated there are denoted by a subscript s. We allow for the possibility that 
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Fig. 2.— SED for the flared blackbody disk. At mid-IR wavelengths, Ly oc At longer 

wavelengths, (xv^. 




the gas temperature in this layer, Tgg, may be smaller than the dust temperature, T^s- 
The rest of the disk is lumped together as the interior. Variables evaluated there carry 
a subscript i. We assume that in steady-state the interior gas and dust temperatures are 
equal; Tgi = T^i = Ti. 

A schematic of how the stellar radiation is reprocessed is illustrated in Figure The 
radiation penetrates the disk to an optical depth of order unity along a slant path inclined 
by an angle a to the surface. Dust particles along the slant path are “superheated” to a 
temperature 



( 11 ) 


which is greater than the blackbody temperature Tbb ~ because their 

absorptivity in the visible exceeds their emissivity in the infrared. The normal optical depth 
of the superheated dust layer at visual wavelengths is ry ~ a. The infrared optical depth is 
smaller still; tir ~ aSa- 


The superheated dust radiates equal amounts of IR radiation into the inward and 


outward hemispheres. The outward directed radiation is similar to that of a dilute 
blackbody. Where the disk is opaque to blackbody radiation at 



(12a) 


the inward directed radiation is thermalized at that temperature. The outer boundary 
of this region is denoted by ath- Just outside ath, the interior is optically thin to its own 
radiation but still opaque to radiation from the superheated surface; in these regions the 
interior temperature is determined by thermal balance to be 



(12b) 


At still greater radii, the encased material is transparent to both its own radiation and to 
radiation from the surface; the internal temperature here is given by 
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Fig. 3.— Radiative transfer in the passive disk. Stellar radiation strikes the snrface at 
an angle a and is absorbed within visible optical depth nnity. Dust particles in this first 
absorption layer are superheated to a temperature T^s- About half of the emission from 
the superheated layer emerges as dilute blackbody radiation. The remaining half heats the 
interior to a temperature Tj. 
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The SED for the radiative equilibrium disk is computed from 

pCLo poo (jq- 

Ly = / da a dzB^(T) , (13) 

Jai J—oo dZ 

where measures optical depth from z to cx) along the axis perpendicular to the disk 
midplane. 


2.3.1. Flat Geometry 


Once again we consider the flat disk, but now under conditions of radiative equilibrium. 
The appropriate expression for the effective temperature is given by equation (Q). For our 
hducial flat disk ath ~ 50 AU. Runs of T^s and Ti as functions of a are displayed in Figure 
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The SED for the flat, radiative equilibrium disk as calculated from equation (P!3D is 
displayed in Figure |^. Its appearance is similar to that of the SED for the flat blackbody 
disk. Over most of the IR, it is dominated by radiation from the optically thick interior. 
The surface layer radiates more than the interior shortward of 6 /im; however, there its 
contribution is hidden by that from the central star. Most of the radiation longward of a 
millimeter comes from the outer, optically thin part of the disk. This accounts for the drop 
of the SED below the extrapolation of the n = 4/3 power law. 


2.3.2. Hydrostatic Equilibrium 

Now we investigate the disk model in which both vertical hydrostatic equilibrium and 
radiative transfer are treated in a self-consistent fashion. The flaring geometry is governed 
by equations (|^) through (|^). In the limit where the grazing angle is dominated by the 
flaring term and the disk is opaque to its thermal radiation, the expression for H/a is nearly 
identical to that given by equation (]T(]|). In other, more optically thin regimes, the flaring 
geometry changes slightly as Tj takes successively different forms. 

A plot of Ti for the flared, radiative equilibrium disk is included in Figure ^ 
Approximate htting formulae for T, and H/a in the three regions of the flared, radiative 
equilibrium disk are as follows. For 0.4 < oau < 84, 
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a (R*) 



Fig. 4.— Temperature profiles of the flat and flared radiative equilibrium models. The dust 
temperature Tds of the superheated layer is independent of disk geometry. Expressions for 
Tds and T* are provided in the text. The flat disk is truncated at Oo = 270 AU (to facilitate 
comparison with the flared models), before the third temperature regime is reached. 
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log u (Hz) 

Fig. 5.— SED for the flat radiative equilibrium disk. Emission from the superheated surface 
is hidden by central starlight. For 30 ^ A(/xm) ^ 1000, optically thick emission from the 
disk interior resembles that from the flat blackbody disk (see Figure |^, but is reduced by a 
factor of 2. Emission longward of 1 mm is optically thin. 
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i7/a ~ 0.17 a^u ; 


for 84 < oau < 209, 


and for 209 < oau < 270, 


H/a ^ 0.59 


Ti ^ 21 K 

oau \ 
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(14a) 

(14b) 

(14c) 

(14d) 


(14e) 

(14f) 


The SED for the flared disk truncated at Oq ~ 270 AU, as computed from equation 
(0), is shown in Figure ^.i Radiation emitted by the surface mirrors that of the interior, 
but is shifted towards and dominates at shorter wavelengths. The near power law rises 
of both disk components of the SED in the mid-infrared are readily obtained following 
the procedure used in the derivation of equation (^. The superheated layer obeys the 
scaling relation ~ 8770 ^X 7 ^; (^Tjs ~ CiL^/2 ~ 0.06(z//10i5Hz)“5/7/^^_ Similarly, the 
opaque interior obeys ~ Sna'^aT/ ~ aL^/2 ~ 0.02(z//10i^ Hz)“^/ 5L*. That the former 
exceeds the latter at hxed wavelength reflects the fact that compared to most of the 
contribution to arises from larger radii and hence larger a. 


3. EXTENSIONS AND REFINEMENTS 


^The disk may terminate in a slow thermal wind since i//a « 1 at oau = 270. 
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Fig. 6.— SED for the hydrostatic, radiative equilibrium disk. At mid-IR wavelengths, 
the superheated surface radiates approximately 2-3 times more power than the interior. 
Longward of 300 /rm, n gradually steepens from about 3 to 3 + /3 as the disk becomes 
increasingly optically thin. 
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3.1. The Planetary Atmosphere Analogy 

Planetary atmospheres are generally transparent to visible radiation above their 
nppermost clonds decks.li Visible radiation is absorbed by particles. Collisions maintain 
thermal eqnilibrinm between gas molecules and the particles. Unit optical depth in the 
thermal infrared is reached above the cloud tops. Pressure broadened molecular lines are 
the main source of infrared opacity.^ The total infrared luminosity, the sum of absorbed 
sunlight and internal heat loss, is emitted as approximate black body radiation. Thus the 
energy balance of a planetary atmosphere resembles that of the blackbody disk.!! This 
analogy probably accounts for the casual acceptance of the blackbody disk model. 

However compelling, this analogy is flawed. T Tauri disks differ from planetary 
atmospheres in a crucial parameter, the vertical gravitational acceleration. For disks, 
gdisk ~ with the orbital angular velocity, and 2 ; the distance from the midplane. For 
our standard parameters, gdisk ~ 0.4 zau ^au much smaller than gpianet ~ lO^cms”^ 

for a ^ Qi- Three consequences of the low vertical gravitational acceleration in disks which 
invalidate the analogy with planetary atmospheres are described below. 

Dust grains settle slowly in disks. The gravitational settling time to the midplane. 


t settle 


rpdVL 


10 ^ 


7 /O.l fim^ 


10 ^gem 2 




(16) 


is independent of a. It is longer than the disk lifetime for particles which most efficiently 
absorb stellar radiation. Moreover, the settling may be slowed by turbulent mixing or 
vertical circulation of the gas. 


The gas pressure in disks is so low that dust grains which are directly exposed to stellar 
radiation reradiate almost all of the stellar energy they absorb. Thus the temperature of the 
superheated grains is determined by radiative balance and is independent of the ambient 
gas temperature. The regulation of the gas temperature is a more subtle issue (cf. § p.2|) . 

The low pressure in disks also means that molecular lines suffer negligible pressure 


^Column densities above the cloud tops are similar to that of the minimum mass solar nebula at oau ~ 1) 
E ~ 10^ gcm“^. 

"^Collision induced dipole transitions of molecular hydrogen dominate in the atmospheres of the outer 
planets. 

®The analogy even extends to the association of the planet’s net luminosity with the accretion luminosity 
of the disk. 
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broadening. Thus the gas opacity is concentrated in narrow, Doppler broadened lines. 
Important molecular sources of opacity must have permitted dipole transitions near the 
peak of the blackbody spectrum corresponding to the ambient temperature T. Furthermore, 
these transitions must connect energy levels which do not lie too far above kT. Water is a 
prime candidate for a coolant at temperatures greater than 100 K. It is abundant, and as a 
consequence of being both a hydride and an asymmetric rotor, it possesses a rich rotation 
spectrum shortward of 300 fim. However, our crude estimates suggest that even if much of 
the cosmic abundance of oxygen is tied up in water, water lines would not cover more than 
a small fraction of the infrared spectrum at any temperature. 


3.2. Energy Balance In The Superthermal Dust Layer 


Stellar energy absorbed by dust grains is lost through emission of infrared radiation 
and collisional transfer to gas molecules. The gas molecules, mostly H 2 and He, gain 
energy in collisions with dust grains and lose it in collisions with molecular coolants such 
as H 2 O. The energy added to the internal degrees of freedom of the coolants in collisions 
with gas molecules and by absorption of dust radiation is lost in radiation to space. Thus 
Tds ^ Tgs ^ Tj., where the subscripts j = {ds,gs,x} denote quantities pertaining to the 
dust, the inert gas molecules, and the molecular coolants, respectively. ^ 

Balancing the rates of energy gain and loss per unit surface area of a dust grain yields 




/ R* 


\ a 




T ^g'^gkiTds ^gs) j 


with Vg the sound speed. For rig -C Ud, where the critical number density 


(16) 
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(17) 


gas-CTain collisions have a negligible effect on Tds and equation ( piBf ) reduces to equation 
(Q).B On the other hand, for rig rid, Tds ~ Tgg. The gas density in the superheated dust 
layer, Ugg, is unknown. It could be as large as no, the value at the disk midplane, or as low 
as Umin, the value obtained by assuming uniform mixing of dust and gas up to the tenuous 


®The symbol refers to the excitation temperature of the coolant transition. 
^We evaluate Vg on the right hand side of equation at temperature Tds- 
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outer layers of the disk. Runs of Ud, uq, and Umin as functions of a are displayed in Figure 
1^. It is seen that even hq falls below rid beyond oau ~ 2. 

Balancing the rates at which inert gas molecules gain and lose energy in collisions with 
dust grains and molecular coolants yields 


ndsCTdiTds - Tgs) ^ n^a^iTgs - T ^), (18) 

where n is a number density, and a denotes a collisional cross-section. The compositions 
and mixing ratios of molecular coolants are uncertain. Photodissociation and incorporation 
into dust grains represent loss mechanisms that can reduce abundance ratios far below 
values based on cosmic abundances. Nevertheless, suppose that the mass in coolant 
molecules is comparable to that in dust grains. Then ndsUd/nxCrx ~ 10“^, the ratio of the 
linear size of a gas molecule to a dust grain. Thus we would not be surprised to discover 
that Tgs Tds - Tgs. 

Estimating excitation temperatures for the transitions in molecular coolants is a 
daunting task. We proceed by making a number of drastic approximations. Each transition 
is assumed to give rise to an optically thick, Doppler broadened line having excitation 
temperature T^. For clarity, we begin by considering the pure radiative excitation of these 
levels. It is not difficult to show that the excitation comes mainly from radiation emitted 
by the superheated dust grains as opposed to that from the inner portion of the disk. 

We specialize further to transitions whose frequencies he on the Rayleigh-Jeans part of 
the blackbody spectrum at temperature Tds- Applying the basic concepts of the escape 
probability formalism for radiative transfer to a two level molecule, we obtain 


Tx ^ aSsTds 
k \aSskTdsJ 


for hu ^ kTx , 

(19a) 

for hu ^ kTx . 

(19b) 


Pure radiative excitation sets a floor on the value of Tx- 

Next we consider the effects of collisions while neglecting radiative excitation. Here 
we suppose that lines of brightness temperature Tx cover a fraction / of the blackbody 
spectrum at temperature Tx- We equate the emission line flux to the rate (per unit disk 
surface area) at which the gas gains energy in collisions with the superheated dust grains: 



18 


a (R,) 

100 1000 10 ^ 



a (AU) 

Fig. 7.— Comparison of possible gas densities in the superheated layer to the critical density 
Uci above which gas-grain collisions might lower the dust temperature. The actual gas density 
could be as high as no, the value at the disk midplane, or as low as n^m, the value for a disk 
in which gas and dust are uniformly mixed everywhere. Except in the innermost 2AU, it is 
evident that the dust likely remains superheated at the temperature given equation (pT|). 
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faT^ ^ ngsVgka{Tds - Tg^) 
For Hgs -C nc 2 , where the second critical gas density 
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nc2 ~-^ci, 

0^6^ o 


we hnd 




12 ^ 2 / 


ds 


In the high density limit Ug S> rid, ~ Tgg ^ T^s ~ Tg with 


a 


1/4 , p X 1/2 


y -|- GEg j \ a 


(-) T. 


( 20 ) 


( 21 ) 


(22a) 


(22b) 


If / -C GEg, equation ( |22b|) reduces to (|^, the appropriate expression for the dust 
temperature in the absence of cooling by molecular line emission. On the other hand, if 
/ S> GEg, (p2b|) collapses to the expression appropriate for a blackbody disk as given by 
equation (|I|), magnified by a factor of about 1//^/'^. 


3.3. Location Of The Superthermal Layer 


We have been assuming that the dust and gas are uniformly mixed throughout the 
disk. Consequently, for our model the gas density is much smaller in the superthermal dust 
layer than in the disk midplane (cf. Fig. ^). Although the timescale for dust grains to 
settle to the disk midplane may be long compared to the disk lifetime [cf. equation (pA]) ], 
the timescale for grains to fall through the superthermal layer (in the absence of advection 
by the gas) is considerably shorter. The latter is given by 


^settle,. 


aVKvPd.k^'^ 


2 


2 /400 cm g 

~ «AU 


r 


K,V 


yr, 


(23) 


where kv is the dust opacity in the surface layer at visual wavelengths. The timescale 
for dust collisions/coagulation to occur in the superthermal layer is l/a times larger than 
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tsettle,s- Thus, the superthermal dust layer is likely to be found closer to the disk midplane 
and to have a lower Ky than our simple assumption of uniform dust-gas mixing implies. 

How might we learn about the level at which the superthermal dust layer sits? One 
possibility is through the observation of emission lines that would form in the gas associated 
with this layer if the gas temperature is comparable to the dust temperature. Perhaps 
most revealing would be observations of the rotation lines of H 2 . Because these arise from 
quadrupole transitions and are therefore weak, they would only be seen if the dust layer 
were located where the gas density is much larger than the minimum value given by the 
assumption of uniform mixing of dust and gas. 


3.4. Dust Albedo 

A nonzero dust grain albedo, A^, at visible wavelengths would reduce the absorbed 
stellar flux by a factor of order (1 — This would lower the disk thickness, but by a 

less impressive factor of order (1 — . 

How small might (1 — be? Values of Ad ~ 0.4 — 0.5 are typical of interstellar 

grains (Draine & Lee 1984). Visual albedos of pure ice particles are much closer to 
unity. And ices are likely to comprise a signihcant fraction of the mass in grains in the 
outer disk. Starting inside oau = 10^, the partial pressures of H 2 O, then NH 3 , and hnally 
CH 4 are expected to exceed their vapor pressures. However, taking the outer planet 
atmospheres as a guide, it seems unlikely that any ice so formed would be pure enough to 
yield (1 - A,)V 2 ^ 0.3. 


3.5. Accretional Heating 

A star which derives its radiant energy from accretion has T* RlaY^^, 

where M* is the mass accretion rate.i The stellar flux impinging on the circumstellar 
disk, aT^{R^,/aYct, is to be compared with the local viscous dissipation rate per unit 
area, (jT^{R^/ aY. Thus in the extreme case that accretion accounts for the entire stellar 
luminosity, viscous dissipation is only competitive with stellar heating for a ^ atr- 

The bulk of the viscous dissipation probably takes place deep inside the disk. Provided 
the disk is sufficiently opaque, this could result in midplane temperatures being higher 


^Provided the accretional energy is radiated over the entire stellar surface. 
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than by a factor of order where r is the effective vertical optical 

depth. A concomitant increase in the disk thickness by a factor of changes the snrface 
geometry and may canse some portions of the disk to be shadowed from direct exposnre 
to starlight. Snch effects are likely to be important only for the innermost regions of the 
disk where radiation temperatnres, snrface densities, and conseqnently optical depths are 
highest. For example, aronnd an accretion powered star, the midplane temperatnre of onr 
standard disk is abont 2000 K at oau = 1-5. This yields abont a factor 2 increase in disk 
thickness relative to that of the passive disk. For smaller radii, the enhanced thickness 
wonld be bnffered somewhat by the redaction of opacity associated with the vaporization 
of silicates. 


4. DISCUSSION 

4.1. Comparison with Observations 

4-1.1. Continuum SEDs 

In Fignre |^, we demonstrate that the nearly flat infrared excess of the classical T 
Tanri star GM Anr can be natnrally interpreted as arising from a passive reprocessing 
disk. Parameters for the central star are obtained from Beckwith et ah (1990), and a 
distance of 140 pc to the Tan-Anr clond is assnmed. Approximate best £t parameters for 
the assnmed face-on disk are the same as those of onr hdncial model, with the following 
alterations: f3 = 1.4, = 0.3/rm, Sq = 3 x lO^gcm”^, a,, = 4 x lO^i?* = 390 AU, and 

= 680i?, = 5.7 AU. 

The inner cnt-off radins is chosen so that the model matches the observationally 
determined SED at 12 /rm; a smaller cnt-off radins results in excessive emission. However, 
it seems likely that the central hole is an artifact of htting a face-on model to observations 
of a disk viewed at a non-zero inclination angle, i. Millimeter-wave aperture synthesis 
images in CO and Hubble Space Telescope (HST) images in scattered stellar light reveal a 
flared circumstellar disk in Keplerian rotation (Koerner 1997), and suggest that i ~ 60°. 
It is possible that the flared outer “wall” extincts emission from the inner disk, diluting 
the 12 fim flux and mimicking the signature of an inner hole. We reserve a more detailed 
discussion of the effects of non-zero inclination angle to a future study. 
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A (miti) 



Fig. 8.— Observed SED of GM Aur (solid points) and accompanying hydrostatic, radiative 
equilibrium (face-on) disk model. Fit parameters are as follows: S = 3 x g cm“^, 

f3 = 1.4, Tp = 0.3 pm, tto = 390AU, and a* = 5.7AU. The derived inner cut-off radius is 
likely an artifact of htting a face-on model to a disk observed at substantial inclination. See 
the text for a discussion. 
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4-1.2. Spatially Resolved Broadband Observations 

Pioneering observations which marginally resolve disks at mm wavelengths are reported 
in Lay et ah (1994) and Mundy et ah (1996). Fits of elliptical Gaussians to the 
A = 0.87mm brightness distributions of HL Tan and L1551 IRS 5 yield semi-major radii at 
half-maximum brightness of 60 and 80 AU, respectively (Lay et ah 1994). Similar results 
were obtained at A = 2.7 mm for HL Tan, with the semi-minor radius marginally resolved 
to be ~ 30 AU (Mundy et ah 1996). 

Interferometric observations of the nearest T Tauri disks at A ~ 10 pm with 
10 milliarcsecond resolution might separate the contributions to the SED from the 
superthermal dust layer and the disk interior. Running integrals for and against 
a are displayed in Figure p. Most of the 10 pm emission inwards of 1 AU originates from 
the disk interior, whereas the contribution from the superthermal layer is localized in an 
annulus centered on ~ 10 AU. Similar results hold for the 20 pm flux. At wavelengths 
longer than ~ 100 pm, the disk interior dominates the emission. 


4-1.3. Spectral Dust Features 

Dust grains possess spectral resonances. Among the best-known are those at A = 9.7 pm 
and 18 pm which are thought to arise from silicates (see, e.g., Draine 1995); the former 
has been observed in both emission and absorption in the spectra of T Tauri stars. For 
example, Cohen & Witteborn (1985) identihed 10 pm emission features in 24 stars and 
10 pm absorption features in another 7 in a spectrophotometric survey of 32 T Tauri stars. 
The emission features evince line to continuum ratios of about 1.2-3.0 to 1. 

Emission features are a natural consequence of the superheated dust layer for disks 
viewed nearly face-on. Line to continuum ratios should be of the same order as the 
percentage increase in dust emissivity associated with the spectral feature. To illustrate the 
effects of spectral resonances, we have employed an emissivity prohle motivated by the data 
for “outer-cloud dust” as given by Mathis (1990). The resulting (face-on) SED displayed 
in Figure 0 exhibits emission features at 10 and 20 pm. 

Absorption features may be associated with disks viewed nearly edge-on. Support 
for this hypothesis comes from two observational hndings. These features appear 
most prominently among the “extreme” flat spectrum sources for which L^, at infrared 
wavelengths matches or exceeds Liy in the visible. There is a positive correlation between 
the strengths of these features and the fractional linear polarization at optical wavelengths. 
Both hndings are plausibly the consequence of the extinction of optical radiation in a hared 
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Fig. 9.— Running contributions to the SED from the superthermal dust layer (S) and 
the disk interior (I) at wavelengths for which the disk may be increasingly resolved by 
interferometers. For a given disk component, Li.(< a) is given by equation (|T^, with Oo 
replaced by a. At A = 1 and 3 mm, the contributions from the surface layer are negligible 
and not shown. 


CO 


o 

X 

I— 

o c 

1 

n ^ 1 ' 1 ^ ^ 

10 /i,m 

S 

n ' 1 ' 1 ' 1 

20 /i,m 

S / 

5x10"^ 

I - 

I / 

/ 

0 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

' 

1 1 1 1 1 1 1 

1 1 1 1 1 1 1 

5x10"® 

III 

1 7 

1 mm ! 

I / ■ 

III 


3 mm j 


1 /■ 

10"® 

/ 

— / - 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


/ 

/ 

5x10"'' 

/ 

_ / _ 

/ 

/ 

/ 

/ 

/ 

/ 


/ 

/ 

/ 

/ 

/ 

0 

J_^_1_1_1_^^_ 

J_1_1_1_1_1_1_ 










25 


disk viewed at large inclination angle. 


4-1 ■4- Temporal Behavior 

Temporal variations would provide another diagnostic of circumstellar disks. Of 
relevance here is the claim by Moriarty-Schieven & Butner (1997) that the submillimeter 
and millimeter fluxes from the T Tauri binary GG Tan increased by factors of order 2 
between 1992 and 1994. The cause of this “radio-wave flare” has not been identihed. A 
plausible hypothesis is that it resulted from enhanced disk heating associated with a burst 
in luminosity originating near one or both components of the central binary. This leads us 
to consider relevant timescales for the radiative and hydrostatic response of the disk. 

Seven different timescales come into play. They are: (1) the timescale over which 
superheated dust grains in the surface layer equilibrate with the ambient stellar radiation 
held, 


tds ~ 


rpdkTdsO^ 

PdRIcrT^ 


~ 0 . 0201(5 s, 


(24) 


where pd ~ 10/ig is the mean molecular weight per degree of freedom in a dust grain; (2) 
the light-travel timescale from star to disk; 


to ~ ~ 5 X 10 oau s ; 

c 


(3) the photon diffusion timescale through the disk,E 


^diff 


EiKvTjh 2 X 10^ 


9/14 

®AU 


(4) the thermal timescale for the dust if it is decoupled from the gas,0 


(25) 


(26) 


Srffc / 1 \ 2x 10^ J (1 + 10-"ai(('") foroAu^Sd 

Pd(TTf + EiKvT.) ^ alS" ^ \ (1 + (^)"''') for 84 ^ oau ^ 270 


( 27 ) 


®This timescale is only relevant where the disk interior is opaque to its thermal radiation. 
^°In evaluating the timescale at the outermost radii, we have set Ti = 21K. 
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Fig. 10.— SED for the hydrostatic, radiative equilibrium disk using a grain emissivity 
prohle motivated by data from Mathis (1990). For wavelengths shorter than 0.3/rm, our 
assumed emissivity is unity; longward of 0.3 yum, it obeys a (single) power law relation 
Sa = (0.3 yum/A)^'"^, on which are superposed two Gaussians centered on 10 and 20y(im, 
having amplitudes that are 3 times their local continuum emissivity and FWHM equal to 3 
and 9 yum, respectively. 
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where 10 is the surface density of dust; (5) the thermal timescale for the gas/*^ 


Efc / 1 \ 2 X 10® / + 10-^aS^^) for uau ^ 84 

V ^ SiKv^) ^ \ (1 + for 84 ^ uau ^ 270 


(6) the timescale for the dust temperature to relax to the gas temperature, 


trelax 


rpd fJ-grpd 

- r\j - 


10 ^aius; 


(29) 


(7) the dynamical timescale over which the disk adjusts to departures from hydrostatic 
equilibrium, 


tdyn ^ lAaXsJi . (30) 

How rapidly might the SED vary in response to changes in the luminosity of the central 
star? Since tds ^ contributions from the surface layer are limited by t^. Those from the 
interior are limited by tgi, since Ueiax ^ tdi and tdi/f tgi^ To relate these response times 
to timescales for variation at a hxed wavelength. A, consnlt Figure 


4.2. Unresolved Issues 

Our investigation leaves many nnresolved issues. 

Is a disk in radiative and hydrostatic eqnilibrinm dynamically stable? 

How does the SED depend npon disk inclination? 

How mnch of the thermally-emitted spectrnm is covered by molecular lines? Do the 
lines appear in absorption or emission? Which molecules are they associated with? 

How high does the layer of superthermal dust grains sit above the disk midplane? 
What is the gas temperature in this layer? Is gravitational settling of dust grains seriously 
impeded by turbulent mixing or by vertical circulation? Or does radiation pressnre from 
the central star dominate all the forces, so that the superthermal grains are driven into 
denser environments and the disk geometry flattens? 


^^The former inequality is only marginally satisfied at ao- 








What is the albedo of the circumstellar dust? How does scattered stellar radiation 
compare to thermal emission in the near IR? 

Does active accretion signihcantly affect the SEDs, either by preventing the settling of 
dust grains or by thickening the inner regions of the disk? 


We thank Anneila Sargent for providing the data on GM Aur and other T Tauri stars, 
Shri Kulkarni for suggesting that temporal variability of disk emission can provide a further 
means of probing disk characteristics, Dave Koerner and Karl Stapelfeldt for informative 
discussions regarding GM Aur, and Steve Beckwith for a thorough and thonghtful referee’s 
report. Financial snpport for this research was provided by NSF Grant 94-14232. E. G. 
acknowledges snpport from an NSF Gradnate Fellowship. 
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